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In the last few years it has become possible to study low energy tttt scattering 
in kaon decays to three pions, thanks to the high statistics measurement of 
K ± —> tt^tt^tt® decays performed by the NA48/2 experiment at the CERN 
<D \ SPS. At the tt + tt~ threshold, the tt°tt° mass spectrum exhibits a Wigner-cusp, 

from which the 5-wave tttt scattering lengths are extracted with high preci- 
sion. This measurement is complementary to the extraction of the scattering 
lengths from K e ^ decays, which is also performed by the NA48/2 experiment. 

^- : 1 Introduction 

q ■ In the low-energy regime, the perturbative description of the strong interac- 

! tion breaks down, as the strong coupling constant becomes of 0{1). Chiral 

Perturbation Theory (ChPT) is an effective theory, which circumvents this 
problem by making use of the chiral symmetry of the theory in the limit of 
vanishing quark masses. Spontaneous breaking of the chiral symmetry gener- 
ates 8 pseudo-scalar Goldstone bosons, among them pions and kaons. They 
obtain their small but non-zero masses by the additional symmetry break- 
ing of non- vanishing quark masses. In the framework of ChPT, the values 
of the iso-spin and 2 5-wave tttt scattering lengths a and a 2 are directly 
connected with the size of the chiral condensate and the pion mass [1]. The 
scattering lengths can accurately be predicted to aom 7r + = 0.220 ±0.005 and 
a 2 m 7r + = — 0.044±0.010 [2]. Thus, precise measurements of the tttt scattering 
lengths are a crucial test of ChPT. 

Previous measurements have traditionally been performed in the semilep- 
tonic decay K ± — ► tt + tt~ e ± is(V) (K e ±). An early measurement by the Geneva- 
Saclay experiment analyzed 30000 events [3] . More recently, the BNL exper- 
iment E865 has measured (a — a 2 )m 7r + = 0.258 ± 0.013 from about 400000 
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Figure 1: Set-up of the NA48/2 experiment. 



K e 4 events [4]. Another recent determination of the scattering lengths has 
been carried out by the DIRAC experiment at CERN from the lifetime of 
pionium atoms. They obtain (a — a 2 )m 7r + = 0.264 from an analysis of 
a part of their data [5] . 

Here, new measurements of the NA48/2 Collaboration are reported. In 
K ± — > 7t ± tt 7t decays, a Wigner-cusp from tttt rescattering in the decay 
amplitude was discovered at m(7r°7r°) = 2m^. This allows a very precise 
determination of a and a 2 with a completely new method. In addition, the 
NA48/2 Collaboration has also determined a and a 2 in K e4 decays, with a 
greatly improved precision with respect to previous measurements. 

The NA48/2 experiment has been taking data in the years 2003 and 
2004. From a 400 GeV/c proton beam, positive and negative kaons with a 
momentum bite of px = (60 ± 3) GeV/c were simultaneously selected by 
a system of achromats (see Fig. []]). A detailed detector description can be 
found in [6]. The main aim of the experiment was the search for direct CP 
violation in decays of charged kaons into three pions. The trigger therefore 
was designed to efficiently select events with three charged tracks as well as 
K ± — > 7T ± 7r°7r events. In total, about 2 x 10 9 three-track events and about 
90 x 10 6 K ± —> 71^71^71® events were recorded in the two years of data-taking. 
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Figure 2: Distribution of m£ 07r o of the NA48/2 K± -> 7r ± 7r°7r° data. Left: 
Full kinematic region. Right: Zoom into the cusp region. The arrow marks 
the (2m 7r +) 2 threshold. 



2 Wigner-Cusp and tttt Rescattering in 

K ± — > 7r ± 7r°7r° Decays 

With the unprecedented statistics collected in the channel K ± —> ^Tr ^ , 
it was possible for NA48/2 to precisely measure the distribution of the 7r°7r° 
invariant mass. The original aim of this analysis was the search for the 
decay K ± — ► 7r ± (7r + 7r _ ) at o m with the subsequent decay of the pionium atom 
(7r + 7r - ) a tom into 7r°7r°. It came as a surprise, that the 7r°7r° invariant mass 
distribution exhibited a distinct and clearly visible cusp at the tt + tt~ mass 
threshold (see Figure [2). 

After its discovery, the effect was explained by N. Cabibbo [7] as Wigner 
cusp, arising from 7r + 7r~ — > 7r°7r° rescattering from the decay K ± — > 7r ± 7r + 7r _ Q. 
As consequence, a rescattering amplitude J\4i from K ± — > has to 

be added to the unperturbed K ± — > ^n^n® amplitude A^o- Below the cusp 
point ttwo = 2m 1T + ) A4i is real and negative, thus explaining the observed 
deficit of events. In the one-loop approximation, Aii turns imaginary above 
the cusp. In this region, the effect on m^o^o is much smaller and less visible. 
The cusp strength is proportional to (a — a 2 )m 7F +, the effect therefore allows 
the extraction of the tttt scattering lengths. 

1 Independently, the effect had already been predicted in tttt scattering much earlier by 
U.-G. MeiBner and collaborators [8]. 
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Figure 3: Diagrams of two-loop 0(a?) corrections (from ref. [9]). 



Second order calculations were performed by Cabibbo and Isidori [9]. 
These include two-loop level 0(a?) corrections (Fig. [3]), which allow the sep- 
arate determination of the scattering length a 2 . Also other rescattering pro- 
cesses as 7r°7r° — ► 7r°7r°, 7r + 7r° — > 7r + 7r°, etc. are considered, and all K — > 3n 
decays as e.g. K L — > 7r 7r°7r° are covered. From this, the theoretical error on 
the extraction of a — a 2 is about 5%. To reach a higher level of precision, 
0(af) and radiative corrections have to be taken into account. 

Using the Cabibbo-Isidori model, the NA48/2 collaboration has per- 
formed a fit to the distribution seen in —> 7T ± 7r°7r [10], following 
a previous analysis on a fraction of the data [11]. Free fit parameters were 
the scattering lengths a — a 2 , a 2 , the linear and quadratic slopes g and h 
of the Dalitz plot u distribution, and the normalisation. The region around 
the cusp point at rn% = (2m 7r +) 2 was excluded from the fit to avoid an 
influence of possible pionium formation. The result of the fit is 

(oo - a 2 ) = 0.261 ± 0.006 stat ± 0.003 syst ± 0.001 ext ± 0.013 theory 

a 2 m^ = -0.037 ±0.013 stat ±0.009 syst ±0.002 ext . (1) 

The result is still preliminary. The residuals of the fit are shown in Fig. [H 
When the correlation of p = —0.92 between both values are taken into ac- 
count, the result on the scattering length a is 

a m, + = 0.224 ± 0.008 stat ± 0.006 syst ± 0.003 ext ± 0.013 theory . (2) 

The main contributions to the systematic uncertainty come from the shower 
simulation in the calorimeter, the trigger efficiency, and the dependence of 
the Dalitz plot v parameter. The external error comprises the ratio between 
the K + — > 7r + 7r + 7r~ and K + — > 7r + 7r°7r° decay widths [13]. The theoretical 
uncertainty is estimated due to the neglected 0{af) and radiative corrections, 
as mentioned above. 

Using analyticity and chiral symmetry, a and a 2 are directly related with 
each other [14]. With this constraint, the experimental error is reduced sig- 
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Figure 4: Pulls of the fit to the rn% distribution. The region between the 
vertical lines was excluded from the fit. 



nificantly and the NA48/2 result becomes (a — a 2 ) = 0.263 ±0.003 sta t ± 
0.001 syst ± 0.001 ext ± 0.013 theory . 

When including the region around the cusp into the fit and allowing 
for an additional pionium component in the fit, NA48/2 observed an excess 
of (1.8 ± 0.2) x 10 -5 events per K ± — > ^n+n - decay at the cusp point 
(see Fig. [4]). This is a factor of 2.3 ± 0.3 larger than the predicted rate 
of 0.8 x 10 -5 for pionium formation in K ± — > ^tt+tt - [15]. However, it 
has recently pointed out [16], that electromagnetic corrections to final state 
interactions have to be taken into account, which might explain the observed 
excess over the prediction. 

3 Rescattering in Kl — > 7r 7r°7r° decays 

In K L — > 7T°7T°7T events, a similar cusp at m^o^o = is expected [12]. By 

computing the matrix elements for K + / ' K L — > 3tt at the cusp point from the 
measured partial widths and Dalitz plot slopes [13], the visibility is expected 
to be about 13 times smaller in K L — ► 7r 7r°7r° compared to K ± — > ^n ^. 

The NA48 collaboration has analyzed about 88 million —> 7r 7r°7r° 
events, recorded in a special data-taking period in 2000 with no material 
in front of the electro-magnetic calorimeter. The data are practically back- 
ground-free. The distribution of squared 7r°7r° invariant mass is shown in 
Fig. [5] (left). It should be noted, that in the region of low invariant mass 
around the 7r + 7r~ threshold the plot contains at maximum one entry per 
event, and therefore no double-counting occurs. When dividing the data with 
events from a Monte Carlo simulation without tttt rescattering, clear evidence 
for a change of slope near the cusp point is observed (Fig. [5] (right)). 
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Figure 5: Distribution of ra^o of the NA48 K L 7r 7r°7r° data. Left: Full 
region. Right: Data divided by Monte Carlo simulated events without rescat- 
tering (arbitrary normalisation). The arrow marks the (2m 7r +) 2 threshold. 



4 Rescattering in — > t^tt+tt decays 

Effects from tttt rescattering should also occur in decays to charged pions 
as in K ± — > ^n+n' . In this case, 7r + 7r~ — > 7r + 7r~ and 71^71^ — > 71^71^ 
contribute. Here no cusp is possible inside the physical region, but the Dalitz 
plot distribution is modified. The K —> 3tt matrix element, as a function of 
the Dalitz plot variables u = (s 3 — <s )/m 2 + and v = (si — 5 2 )/m^ + , is given 
by [13] 

|.M| 2 = l + ^^ + /i^ 2 + fc^ 2 +••• , (3) 

with the slope parameters g,h, and k. 

NA48/2 has performed a combined fit to their K ± — > tp^^ and X ± — > 
7r^7r^7r — data, using the non-relativistic effective field theory of the Bern- 
Bonn group [17]. Free parameters in the fit were the normalisations, — a 2 , 
and the slope parameters g ) h, and k of K ± —> ^tt+tt - and g and h of 
— > ^tt ^ . The fit had a nearly perfect x 2 of 757.1/757 d.o.f. and 
resulted in significantly shifted values of the slope parameters g, /i, and k 
with respect to a fit with no rescattering taken into account. However, in the 
latter fit the x 2 value was equally good (516.9/517). This means, that on the 
one hand nn rescattering has a significant impact on the values of the Dalitz 
plot slopes. On the other hand, there is practically no possibility to extract 
the scattering lengths from the K ± — > ^tt^^ Dalitz plot slopes alone. 
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Figure 6: Distributions of the five Cabibbo-Maksymowicz variables in the 
NA48/2 2003 K e4 data. 



5 Measurement of tttt Scattering Lengths in 
K e 4 Decays 

An independent and complementary method for the determination of tttt 
scattering lengths is the analysis of K ± — > tt + tt~ e ± v(v) {K e ^) decays. The 
K e 4 decay is a rare decay with a branching fraction of about 4 x 10 -5 [13]. 
Its amplitude depends on the two complex phases 5 and Si, which are the 
S and P wave tttt scattering phase shifts for isospin 1 = 0. In K e ^ decays, 
their difference 8 = S — Si can be measured as a function of m^. 

In the 2003 data-taking period, the NA48/2 collaboration has collected 
677510 K e 4 decays with very small background contamination of only about 
0.5% [18]. An analysis was carried out in the five independent Cabibbo- 
Maksymowicz variables: these are the squared invariant dipion and dilepton 
masses and s e , the angles 9 n and 9 e of tt + and e + with respect to the tttt 
and ev directions in the tttt and ev rest frames, respectively, and the angle 
(j) between the tttt and ev decay planes. The distributions of these variables 
are shown in Fig. [6l 

NA48/2 has performed a combined fit to the decay form factors and 
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Figure 7: Measurements of the phase difference 5 S — 5 P as function of 
K e 4 decays. The line shows the best fit to a and a 2 of NA48/2. 



m 



the phase shift difference 5 as a function of m^. The results are shown in 
Fig. H together with the earlier measurements of the Geneva-Saclay [3] and 
BNL E865 [4] experiments. The results are in good agreement with each 
other (except perhaps the highest data point of E865). From the phase shift 
measurements, the tttt scattering lengths can be extracted. Using dispersion 
relations (Roy equations) and data above 800 MeV, a 2 is related to a (the 
so-called Universal Band). From a one-parameter fit, NA48/2 obtained 

+0.018 

a m, + = 0.256 ± 0.006 stat ± 0.002 syst , (4) 

— U.U1 I theo 

where the theoretical error corresponds to the width of the universal band, 
which is given by the experimental uncertainties of the measurements above 
800 MeV. This result implies a 2 = -0.031 ± 0.001 exp ± ^012 theo' When 
both a and a 2 are left free in the fit, the results are 

a m, + = 0.233 ±0.016 stat ±0.007 syst , 

a 2 m^ = -0.047 ±0.011 stat ±0.004 syst . (5) 

It has recently been realized, that effects from isospin symmetry breaking 
(rrtd ^ m u ) are not negligible and do have a significant impact on the value 
of the phase shift difference 5. A preliminary ChPT calculation gives a shift 
of S between 12 and 15 mrad for > 300 MeV/c 2 [19]. Taking into 

account these corrections, the value of a shifts by ~ —0.02. Note, that 
these corrections also have to be applied on the data of the previous K e ^ 
measurements. 
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no isospin breaking corrections with isospin breaking corrections 




Figure 8: Comparison of results on the scattering lengths a and a 2 from 
the NA48/2 and DIRAC experiments without (left) and with (right) isospin 
breaking corrections applied to K e ^. Shown are also the universal band (black 
band), predicted using the Roy equations and higher energy data, and the 
ChPT prediction (narrow blue band). 



6 Conclusions 

In the past year, there has been significant progress in the extraction of 
tttt scattering lengths from kaon decays. The NA48/2 experiment has dis- 
covered a Wigner-cusp in the m^o^o distribution of K ± — > tt ± tt°tt° decays, 
arising from tttt rescattering in this channel and thus allowing the precise 
extraction of the scattering lengths ag and a 2 - At the moment, the uncer- 
tainty is dominated by the neglect of radiative corrections and higher order 
terms. Given the on-going theoretical work, this is expected to improve in the 
near future. In addition, evidence for tttt rescattering in K L —> tt°tt°tt events 
from K L —> tt+tt~tt° has been found, using NA48 data from a K L run period. 
The more traditional method for determination of ag and a 2 is the angular 
analysis of K e ^ decays, where the NA48/2 collaboration has presented a new 
analysis with much larger precision than previous experiments. 

The constraints on a and a 2 from the NA48/2 measurements, together 
with the recent a — a 2 determination from pionium lifetime by DIRAC [5] , 
are shown in Fig. [8l When taking into account the corrections for isospin 
symmetry breaking in K e ^ there is very good agreement within the results 
and with the ChPT prediction. Further improvements are expected in the 
near future from the analysis of the complete NA48/2 K e ^ and DIRAC data 
sets, as well as from a better theoretical understanding on the extraction of 
scattering lengths from K ± — > ^tt ^ and K e ^. 
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